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ABSTRACT 

The fraction of matter that is in the form of baryons or dark matter could have spatial fluctuations in 
the form of baryon-dark matter isocurvature fluctuations. We use big bang nucleosynthesis calculations 
compared with observed light element abundances as well as galaxy cluster gas fractions to constrain 
cosmological variations in the baryon fraction. Light element abundances constrain spatial variations 
to be less than 26-27%, while a sample of "relaxed" galaxy clusters shows spatial variations in gas 
fractions less than 8%. Larger spatial variations could cause differential screening of the primary 
cosmic microwave background anisotropics, leading to asymmetries in the fluctuations and ease some 
tension with the halo-star 7 Li abundance. Fluctuations within our allowed bounds can lead to "B- 
mode" CMB polarization anisotropics at a non-negligible level. 
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1. INTRODUCTION 

Is the cosmological fraction of matter that is in the 
form of baryons a universal constant? There is no com- 
pelling theoretical model that would predict it to be ob- 
servably non-constant (although such an effect would be 
fairly easy to "explain"), but there is little empirical ev- 
idence that it is indeed a constant to high precision. 

The cosmological gravitational potential fluctuations 
are observed to be at the level of 10~ 5 , but this measures 
fluctuations in the sum of baryonic and non-baryonic 
mass. The fraction of matter in the form of baryons could 
have much larger long wavelength fluctuations with little 
observable imprint in the cosmic microwave background 
(CMB). Such a modulation is a baryon-cdm isocurvature 
mode, and its imprint on the CMB is a second order ef- 
fect. On scales larger than the sound horizon at recombi- 
nation there is almost no effect, while on scales that are 
sub-horizon at last scattering there would be an addi- 
tional modulation of the acoustic peaks due to a varying 
sound speed. The extra pressure due to the baryons is 
negligible compared to the photon pressure. 

In this paper we investigate other effects of the varia- 
tion of the baryon fraction, and we find that large scale 
modulation is not only allowed but would actually alle- 
viate some current mild tensions in the standard cosmol- 
ogy- 

In particular, we below investigate constraints from 
big-bang nucleosynthesis (BBN) in light of observed light 
element abundances, observations of galaxy cluster gas 
fractions, and modulation of the optical depth to Thom- 
son scattering of CMB photons. We also investigate the 
expected contribution to CMB polarization anisotropy 
and discuss possible other probes. 

2. BIG BANG NUCLEOSYNTHESIS 

The light-nuclide yields from BBN are functions only 
of the baryon density ps of the universe under the stan- 
dard assumptions of uniform entropy per baryon, only 
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standard-model particles, neutrino-antineutrino asym- 
metries not enormously larger than the corresponding 
baryon asymmetry, and no late additions of entropy 
([Schramm fc TurnerllI998l ). Observations of extragalac- 
tic deuterium provide a particularly tight constraint on 
Pb based on the steep dependence of the deuterium yield 
on this parameter. Using VIbH 2 = 8itGpb/3 to ex- 
press the constraint in terms of the fraction fig of the 
closure density provided by baryons, the average extra- 
galactic D/H (= ratio by number of deuterium to hydro - 
gen) gives n B h 2 = 0.0213 ± 0.0010 (jPettini et al.ll2008ft . 
(H = hx 100 km s _1 Mpc" 1 is the Hubble "constant.") 
CMB anisotropy measurements with the WMAP satel- 
lite n ow imply n B h 2 = 0.02273±0.00062 ijDunklev et all 
2009), in agreement with D/H. In these terms, the baryon 
fraction of the total mass density pu is Jb = Pb/pm- 

There is an extensive literature on the effects 
of inhomogeneities in the bary on distrib u tion d uring 
BBN, going back t o at least IWagonerl (l973) and 
lEpstein fe PetrosiarJ (|1975ft . For a review of baryon- 
iiihomogene ous and other non-standard BBN models up 
to 1993, see lMalanev fe Mathewsl (| 1993ft . 

Baryon-inhomogeneous scenarios of BBN may be di- 
vided into three cases, depending on the typical scale of 
the fluctuations: scales that are comparable to particle- 
diffusion horizons during BBN (such that the physics of 
BBN is affected), larger scales that are still small com- 
pared to the scale of our Galaxy (such that observations 
within our Galaxy are a blend of many different vol- 
umes that individually underwent nearly homogeneous 
BBN), and large scales today (such that cosmological 
scales sample distinct regions of uniform baryon density) . 

In the first case, which dominates the literature on in- 
homogeneous BBN, the inhomogeneities are comparable 
in length scale to the neutron diffusion length during 
BBN (~ 4 x 10~ 5 comoving pc at the start or ~ 0.08 
comov ing pc at the end of BBN; iKurki-Suonio et "all 
(|1997h ). Such short-length-scale inhomogeneities only 
have strong effects if the density contrasts are large (ra- 
tios of ~ 10 6 between high- and low-density domains), 
as was suggested in the 1980s on the basis of a strongly 
first-order QCD phase transition (|Wittenlll984f L A first- 
order phase transition of this kind has since fallen out of 
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favor but is not completely ruled out (|Bovanovskv et all 
2006). The corresponding BBN models tend to over- 
produce lithium and underproduce deuterium relative to 
current me asurements ; for a recent assessment of con- 
straints, see lLara et al.1 (|2006| ). This is not the case that 
we consider in the present work. 

Inhomogeneities longer than the diffusion length re- 
sult in BBN at each location occurring indepen- 
dently at different pr. These were investigated by 
lEpstein fe Petrosianl (|1975f ) . and the most elaborate cal- 
culations were performed with a view toward constrain- 
ing isocurv ature baryon fluctua t ions a s seeds of structure 
formation (iJedamzik fc Fuller! 119951 : ICopi et all 119951 : 
iKurki-Suonio et al.lll997l ). These authors assumed that 
any single observation is a blended sample containing 
matter drawn from the full distribution of pg ; the result- 
ing effect on BBN is a "smearing out" of the standard 
BBN predictions of abundance versus mean baryon den- 
sity. Calculations in the three works last cited mostly 
excluded perturbations exceeding the Jeans mass at re- 
combination (which would not be observable sites of low 
metallicity today), with the effect of diminishing the con- 
tributions of high-ps domains. D/H then becomes much 
larger than observed and the primordial 4 He mass frac- 
tion Yp is reduced relative to the homogeneous case. 

The main constraint on blended inhomogeneities arose 
from flattening of the Li/H versus ps curve, which lim- 
its mixed inhomogeneities to Sps/pB ^ 1- If the Li/H 
observed in low-metallicity halo stars is undepleted, the 
primordial Li/H is constrained to lie near a minimum 
of the Li/H versus pb curve of standard BBN. Any 
mixing with regions of pp far from the minimum in- 
crea ses Li/H unacceptab ly. [See in particular Figs. 1- 

3 of IKurki-Suonio et al.1 (jl997t ) for these results). The 
difference of the WMAP- and D/H- inferred pp from the 
Li/H minimum obviously complicates this inference. 

The type of inhomogeneity to be addressed here is on 
length scales that did not mix during structure forma- 
tion, so that abundance measurements are not blends. 
This possibility was noted in earlier work, particularly 
in light of large differences betwee n the early claims 
of ex t ragalactic D/H m easurements ()Kurki-Suonio et all 
Il997t ICopi et aT] [l998). It was not extensively pursued 
because most observational constraints (all Li/H and 
3 He/H, and all D/H before 1995) were within the Galaxy 
and presumed to constitute a single well-mixed sample, 
and because the higher claimed D/H values were even- 
tually discounted. 

We now examine the implications of the existing light- 
nuclidc abundance data for unmixed inhomogeneities. 
The emphasis in the earlier literature was on how in- 
homogeneities affected the BBN upper bound on flp ; we 
take ttp as a more or less settled quantity and ask how 
large fluctuations away from it can be. 
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The strongest BBN probe of pp is deuterium, with 
D/H oc pg 1 ' . It also constrains site-to-site variations, 
since observations of deuterium believed to be primordial 
span a range of redshifts from about 2.5 to 3.6 at varying 
directions on the sky; the corresponding length scales are 
comparable to the present cosmological horizon. 

In the dashed curves of Fig. [TJ we show likelihood as 
a function of fish 2 derived from seven measurements 




Fig. 1. — (Left) Likelihoods for Slfl/i based on measured abun- 
dances in various locations, shown against a band indicati ng the 
la interval from WMAP5. The BBN yields of lBurles etUI dSjOulT) 
were used, assuming Gaussian distribution of observational errors 
only; error estimates for the BBN calculation are not included. 
(Right) Galaxy cluster gas fraction probability distributions rela- 
tive to mean of sa mple, assuming G aussian probability distribu- 
tions for sample of Allen et al. (2008). 

of D/H in quasar absorption- line systems. These 
are t he results t a bulate d in iKirkman et"aT] (I2003D 



[from iTvtler et alJ (fl9 96) ; 



O'Meara et all 



Buries & Tvtlcr (1998a 
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plus t hose of lO'Meara et al.1 
andlPettini et al.l (120081 ). Following the discus - 
sio n in IKirkman et alP^OOSf ). lO'Meara etail (|2006l) . 
and lPettini et a l. (2008), we have omitted from the figure 
and subsequent analys is the D/H toward Q0347-3819 
(|Levshakov et al.1 l2002f ) and the absorber at z = 3.256 
towards PKS1937-1009 (jCrighton et al.ll2004D (the error 
on the latter being claimed to be underestimated). 

It is well established in the literature that there is 
dispersion within the D/H sample in excess of the re- 
ported uncertainties, so that for example a jackknife esti- 
mate of the error on log D/H is ab out twice that derived 
by the inverse-variance estimate (IKirkman et alJ 120031 : 
lO'Meara et~alll2006l : iPettini et al.ll2008fh It may simply 
be that the errors have been underestimated. Alterna- 
tively, one can accept the published error estimates and 
use the data to limit variation in pb among the locations 
where deuterium is observed. 

2.2. 4 He 

Because the primordial 4 He mass fraction Yp depends 
very weakly on pb, useful constraints require very dif- 
ficult percent-level abundance determinations in sites 
of low metallicity. Operationally, this means observa- 
tions of Hll regions in blue compact galaxies (BCG) o r 
the Magellanic clouds (most recently |Qliye et "ap (|1997h : 
iPeimbert et~aT1 (|2000L[2007| ): Hzotov et al.l (|2007T lL Since 
these observations are extragalactic, they might be sen- 
sitive to large-scale abundance variations. The sam- 
ple of Hn regions well-meas ured for Yp extend s only to 
z ~ 0.05, or - 200ft," 1 Mpc. Ilzotov et al] (|2007f ) recently 
added to the available sample 271 blue compact galax- 
ies identified from the Sloan Digital Sky Survey (SDSS); 
most of these are in the same redshift range as the pre- 
vious sample, but there are tens of objects in the range 
200 to 500 h- 1 Mpc (0.05 < z < 0.13). 

The inferred Yp values are known to suffer from sys- 
tematic difficulties. Disagreement among groups using 
different methods exceeds their claimed errors, and as- 
sessments of the situation are sometimes very pessimistic 
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(|01ive fc Skillmanl [2004) . After a recent revision of the 
atomic data on which Yp rests, some researchers find 
value s in good agreement with WMAP (Pcimbcrt et al. 
l2007fl . while others find Y p too high for a comfortable 
match (jlzotov et al.ll2007f ). 

The systematic problems with Yp should not mat- 
ter for present purposes, because spatial variation of 
Yp should show up as dispersion within any uniformly- 
analyzed data set. However, dispersion is not as simple 
as variation away from the mean Yp, because truly pri- 
mordial Y is not measured anywhere. All Hn regions 
observed for Yp have small but nonzero metal content, 
and the same stars that made the metals also made 4 He. 
Even though systems observed for Yp are chosen for low 
metallicity, inferred values of Yp are almost always based 
on extrapolation to O/H = 0. One either applies a Y 
versus O /H relation calibrated elsewhere or carries out a 
linear regression of Y versus O/H among observed Hn re- 
gions. Despite these complications, no signs of dispersion 
beyond the error bars have been claimed in the modern 
data set. 

In our analysis, we place limits on dispersion within 
the large "H eBCG" samp l e of 9 3 uniformly-analyzed Hn 
regions from llzotov et al] (|2007f ). with Y corrected to Yp 
by subt racting the Y versu s O /H regression slope de- 
rived in llzotov et a l. (2007). We have not included the 
additional SDSS Hn regions from the same publication. 
These were not subjected to the same data reduction as 
the HeBCG sample; they were selected differently; they 
are overall of higher metallicity; their regressed Yp are 
consistent with the HeBCG sample at the 2a level, but 
regressions of the two sets have rather different slopes. 
The errors on the SDSS sample are also large enough that 
a joint regression of bot h sets is still dom inated by the 
smaller HeBCG sample (jlzotov et al.ll2007l ). The con- 
straints on Yp based on each of the 93 Hn regions are 
shown as dotted curves in Fig. [T] 

2.3. 7 Li 

Measurements of primordial lithium arc confined to 
the atmospheres of metal-poor stars in the Galactic halo. 
Since our present interest is in variations of pp on much 
larger length scales, all of these stars together constitute 
a single sample. 

Dispersion of Li/H among the stars of the "Spite 
plateau" - the flat region in the graphs of Li/H versus 
both effective temperature and metallicity - has been ac- 
tively sought ever since the discovery of t he plateau and 
its id entification with the primordial Li/H ([Spite fc Spite! 
1982). This is because star-to-star variation along the 
plateau would be a strong i ndication that the obs erved 
lithium has been depleted (jDelivannis et al.lll993l ). the 
dispersion having arisen from varying amounts of deple- 
tion. Some authors claimed dispersion of 0.04 to 0.1 dex 
(|Thorburnl fl994t iDelivannis et alJ 1 19931 ). but more re- 
cent work has mostly stressed its absence (iRvan et all 
19991: ICharbonnel feJPrimasl 120051: lAsplund et all 120061 : 
Bonifacio et al In the solid curves of Fig. 

U we show the galactic h alo as a single samp le of 
Li/H, using the value s from lAsplund et all (|2006f ) and 
iBonifacio et all (|2007l ). Primordial Li/H values pub- 
lished over a 30-year span all agree within the quoted 
errors. 

These observations at face value indicate that our lone 



sample of primordial Li/H is a factor of three below the 
Li/H predicted by the VLph 2 from WMAP. Spatial varia- 
tions about the WMAP mean would help to explain this 
discrepancy if we live in a region of low pp.- However, 
inspection of Fig. Q] shows that variations in excess of 
50% would be needed for a complete explanation. This 
seems implausible based on constraints derived below. 

2.4. 3 He 

There is no convincingly primordial observed 3 He/H: 
all observations of 3 He are near solar metallic- 
ity, and a strong increase of 3 He/H with time 
is predicted from models of stellar nucleosynthesis. 
Measuremen ts of 3 He/H in variou s locations within 
the G alaxy (iGeiss fe Gloecklerl 11993 iGloeckler fc Geissl 
119961 iBania et al.ll2002t l2007n fail to see any evolution 
of 3 He/H and are all near the WMAP-inferred primor- 
dial value. Presumably some mechanism of 3 He de- 
struction is active, and its outlines seem clea r ()Hoganl 
H995llWasserburg et al.lfl995l : ICharbonnellll995f ). Among 
other things, it causes difficulties for the traditional 
(D+ 3 He)/H constraint applied in the earlier literature 
on inhomogeneous BBN. 

Since all existing 3 He/H measurements are within the 
Galaxy, are near the WMAP-inferred prediction, and suf- 
fer from large uncertainties in the post-BBN evolution, 
they do not provide constraints on variation of pp. 

3. GALAXY CLUSTER GAS FRACTIONS 

Measurements of the galaxy cluster gas fraction can 
provide strong constraints on the cosmic baryon frac- 
tion, provided the galaxy clusters ad equately sample th e 
cosmic baryon fraction distribution (W hite et alJ Tl993). 
Galaxy clusters emit in the X-ray through a combination 
of thermal bremsstrahlung and line emission, allowing 
detailed reconstructions of both density and temperature 
profiles. The electron density can be directly integrated 
to obtain the total gas mass, while the temperature in- 
formation can be added to obtain an estimate of the total 
mass (assuming hydrostatic e quilib rium). 

The sample of I Allen et all (|2008f ) has gas fraction de- 
terminations of 42 galaxy clusters covering the redshift 
range z ~ 0.06 — 1. This is a large homogeneously ana- 
lyzed sample of galaxy clusters that was primarily se- 
lected to include objects that appear to be relatively 
undisturbed ( "relaxed" ) . 

In practice, there are real concerns about the impor- 
tance of the assumption of spherical or biaxial symmetry 
(even triaxial ellipsoidal symmetry is a strong assump- 
tion), deviations from equilibrium, non-thermal pressure 
support, non-equilibration of electron and proton tem- 
peratures, and other issues. However, empirically, this 
method has been found to be relatively robust when ap- 
plied to simulated clusters with many real- world compli- 
cations (|Nagai et al.ll2007T ) . and the restriction to only in- 
clude relatively undisturbed clusters should significantly 
reduce many of these complications. In general, it is 
found that galaxy cluster gas fractions are below the 
cosm ic mean when measu r ed well within th e virial ra- 
dius (jVikhlinin et alJl2006t lAllen et al.ll2008l ). but there 
is little e vidence for c l uster- to-cluster scatter, as seen in 
FigureQ] lAllen et all IpflOSl ) find / gas = 0.1105 ± 0.0016, 
while the WMAP5 baryon fraction is 0.17 ± 0.01. 
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Fig. 2. — Probability distribution of scatter parameter for various 
data sets 

4. CONSTRAINTS ON BARYON FRACTION VARIATIONS 

A maximum likelihood fit was done to each data set, 
varying the mean, m, and intrinsic scatter, s. The fits 
were done to the gas fractions and abundance determi- 
nations. The figure of merit used was 

% 1 

with di the data and at their errors. 

Marginalizing over the mean, an estimate of the likeli- 
hood of the intrinsic scatter was obtained for each data 
set. For the galaxy cluster gas fractions, this was trans- 
lated directly int o scatter in the frac tional baryon frac- 
tion. The fits of iBurles et alj ( |2001l) . ex panded around 
the best fit WMAP5 value for n R /t 2 i|Dunklev et all 
2009), were used to translate the fitted scatter in each 
light-nuclide abundance into scatter in the baryon den- 
sity. Linear fits were applied to the BBN yields in the 
region of interest, providing: 

&n B h 2 ^ AY p 

n B h 2 ~ 0.0087 1 ' 

for the Helium abundance and 

A^fo 2 A[\og 10 (D/H)} 
fl B h 2 0.69 1 ' 

for the Deuterium abundance. 

Results of the likelihood calculations are shown in Fig- 
ure[2] Upper limits on intrinsic scatter in the baryon frac- 
tion were obtained by integrating the likelihoods starting 
at s = until 95% of the probability was enclosed. Con- 
straints from 4 Hc and D/H allow 95% confidence upper 
limits on intrinsic scatter of <7(/b//b) < 0.27 and 0.26, 
respectively. From galaxy cluster gas fractions, this limit 
is only 0.08. (We note that the best-fit s for log 10 D/H 
at s ~ 0.07 (5pb/pb ~ .10) in Fig. [2 ] is in agreement 
with a similar analysis bv lPettini et alT(|2008f ).) 

The scales probed by these measurements vary greatly: 
the 4 He measurements lie at low redshift (less than 0.1), 
the galaxy clusters are at intermediate redshift (less than 
1), and the D/H constraints are at z ~ 2.5 — 3.5. 



The possible selection effects are also very different. 
The 4 He constraints come from a large sample of nearby 
low-metallicity galaxies experiencing star formation to- 
day. The D/H measurements require a strong signature 
of D in a QSO absorption spectrum (hence a large hydro- 
gen column density) and favorable conditions for fitting 
D and H column densities (hence a simple velocity struc- 
ture with minimal nearby absorption and a bright QSO). 

The galaxy cluster sample is selected for being prefer- 
entially "relaxed" in their appearance. This means ap- 
pearing undisturbed and relatively circular. In a universe 
with varying baryon fractions, such selection could be 
problematic: low baryon regions are lik ely to have more 
densi ty fluctuations on small scales (|Eisenstein fc Hul 
1998), and clusters there are likely to form earlier. If 
a sample is selected for clusters that are apparently dy- 
namically more mature ( "relaxed" ) , this sample could be 
skewed toward regions of lower baryon fraction, where 
clusters formed earlier. A sample of galaxy clusters se- 
lected in this way would thus show less scatter than that 
of the universe as a whole, and would be systematically 
biased toward regions of lower baryon fraction. Detailed 
numerical simulations of galaxy cluster formation will be 
required to better understand this important selection 
effect. 

In summary, there is no evidence for cosmological 
variation in the baryon fraction, but the amplitude is 
constrained to be no larger than roughly 10-30%. For 
comparison, baryon-photon or cdm-photon isocurvature 
modes have been constrained to have amplitudes on 
horizon scales at the level of 10~ 5 or smaller (e.g., 
iSollom et alj (|2009f ) and references therein). 

5. DIFFERENTIAL THOMSON SCATTERING 

In recent years, there have been indications of 
a possibl e large-scale asymm e try in CMB fluc- 
tuations (lErikseri et all 12004 IHansen et~aT1 12004 
i Bernui et alJ 120061: lEriksen et all l2007t IHansen et all 
2008: iHoftuft et al J 120091 )7" 

Large scale baryon fraction fluctuations would cause 
such an asymmetry through differential Thomson scat- 
tering. If there are significantly more electrons along 
different lines of sight, this leads to a differential damp- 
ing of primary CMB fluctuations. The observed dipole 
asymmetry in the CMB fluct uations is found to be 
0.072 ± 0.022 [68% confidence, IHoftuft et all (~2009h ] in 
the temperature fluctuation amplitude, and the cur- 
rent best fit for the total Tho mson optical depth is 
0.087± 0.017 (jDunklev et al.ll2009l) . Thus, to explain the 
observed anomaly would require variations of 0.8±0.3 (at 
68% confidence) in the baryon fraction on cosmological 
scales. Such a large variation (but only 2.7c from zero) is 
observed in neither D/H abundances nor galaxy cluster 
gas fractions, but smaller variations could be consistent 
with both the observed asymmetry and the light-element 
constraints. 

Thomson scattering by a region with a fluctuat- 
ing baryon fraction would also lead to polarization 
anisotropies in the CMB. The polarization anisotropy 
can be expressed as a sum o f "E modes" and "B modes" 
(jZaldarriaga fc Seliakl ~997). where E/B refer to curl- 
free and curl components of the polarization field. The 
primary CMB components generated by scalar pertur- 
bations have zero B modes, with the dominant primary 
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contributions to B modes sourced by gravitational radi- 
ation. 

The CMB polarization signature from an inhomoge- 
neous distribution of free electrons during and after 
rcionization has two components: that d ue to differen- 
tial s creening and that due to rescattering (Dvorkin et al. 
2009). The differential screening term arises because 
modulation of the CMB polarization from recombination 
by a direction-dependent e _T ( n ) generates additional B 
and E modes. Although the induced B modes from this 
term are significant at very large angular scales, they are 
subdominant at the intermediate angular scales impor- 
tant for gravitational wave studies. 

The other contribution is due to the scattering of inci- 
dent quadrupole radiation on free electrons, which gen- 
erates both E- and i?-mode polarization if the scatterers 
are inhomogeneous. On small scales (where the spherical 
curvature of the sky can be ignored) the B modes from 
rescattering trace the angul ar fluctuations in the opti- 
cal d epth on the sky CJ T as (|Hdl2000b iDvorkin fe Smith! 
l2009h 

^BB^catt) _ ^.EE^catt) 3 ^ 2 — 2 T e f f /"YTT I a\ 

where Qrms is the RMS quadrupole, which we take to 
be 22 fiK, and T e s is a characteristic optical depth ac- 
counting for the low redshift screening of the generated 
anisotropies. 

Although a primordial baryon-cdm isocurvature mode 
could in principle take any form, with the details set by 
the physics of how it is generated, a natural choice is 
the scale-invariant A: 3 Pb(/c) = constant, where Ps(fc) = 

i 8p B {k) 8p B {k) > integral along the line of sight se- 

N Pb pb ' ° ° ° 

lects a plane of modes in the 3D Fourier space. This 
is simply a slice through the original 3D Fourier space, 
so the statistics of the selec ted modes ar e unchanged 
[the "Limber approximation:" [Kaiser! (|1992l )]. and Pe(fc) 
maintains the same scaling with wavenumber after pro- 
jection. Thus we assume CJ T oc £~ 3 . Figure [3] shows 
the rescattered B-mode polarization due to this effect, 
for total fractional optical depth fluctuations of 20% and 
2%. We show only the result for £ > 50; lower multi- 
poles require a more sophisticated treatment to include 
the effects of the spherical geometry. 

The rescattered B modes in this scenario are well be- 
low the B modes induced from E modes through gravita- 
tional lensing by foreground structure. However, lensing 
estimation techniques could red uce the lensing j3-mod e 
noise by factors of at least 40 ijSeliak fe Hiratal I2004T ). 
Rescattered B modes could consequently be a contam- 
inant in the search for primordial gravitational waves 
using CMB polarization, if the tensor-to-scalar ratio is 
sufficiently low. 

6. FUTURE TESTS 

With current CMB data it should be possible to look 
for variations in the baryon fraction at recombination. 
One could either divide the sky into many segments and 
analyze each section independently ([Hansen et al.ll20C58h 
or use a q uadratic estimator [as was done for lensing 
detection in lSmith et al.l ([2007A ] to map out baryon frac- 
tions on the sphere. 

Large scale structure studies can constrain large scale 
variations in the baryon fraction through baryonic ef- 
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Fig. 3. — The B-modc CMB polarization anisotropy induced by 
scattering with a reionized region containing a power-law baryon 
inhomogeneity with total optical depth variations of 20% (top black 
dot-dashed) and 2% (bottom black dot-dashed), compared with the 
B-mode signal from lensing (top red dotted) and the same curve 
lowered by a factor of 40 as motivated by delensing studies (bot- 
tom red dotted). The induced B modes simply trace the assumed 
optical depth fluctuations with a constant prefactor in the flat-sky 
limit. The red dashed curves show inflation-induced B modes with 
tensor-to-scalar ratios of 0.1 and 0.001. 

fects on the matter transfer function. O ne could use 
galaxy surveys or QSO counts ffhirata 200§) to put limits 
on large scale variations in the matter transfer function. 
This will be complicated by the possibility of the baryon 
fraction fluctuation being correlated (or anti-correlated) 
with the potential fluctuations, but it provides strong 
constraints. There is an additional complication in that 
one needs to understand how the astrophysics depends 
on the baryon fraction. For example, the number den- 
sity of high redshift quasars ([Hiratal |2009| ) is very sensi- 
tive to the matter power spectrum, but it could also be 
extremely sensitive to the baryon fraction. Cooling pro- 
cesses often scale as density squared, so formation times 
could be strongly affected by a varying baryon density. 

Assuming that the varying baryon fraction does not 
aff ect the astrop hysics, we can translate the results 
of iHiratal (|2009l ) into limits on baryon fraction fluc- 
tuations. Baryons aff ect the matter power spectrum 
(|Eisenstein fc Hulll998f ) such that a 1% variation in cr 8 
would be caused by of orde r a 5% c hange in the baryon 
fraction. The constraints of Hirata ( 2009) on large scale 
gradients in the high redshift QSO number density (hori- 
zon scale gradients in In erg less than 0.03 at 99% con- 
fidence) are thus comparable to the galaxy cluster gas 
fraction constraint. 

There is significant room for improvement with a larger 
sample of QSO absorption studies to constrain intrinsic 
scatter in D/H. Deuterium is an excellent measure of 
baryon density; with a modest improvement on current 
samples it is likely that D/H will be the most robust 
constraint on large scale variations in the baryon fraction. 

With future radio interferometers such as the Square 
Kilometer Array, direct measurements of D/H and 
3 He/H will be possibl e over a cosmological vol- 
ume through radio lines (ISigurdson fc Furlanettol 120061 : 
McQuinn fc Switzer![2~009HBagla fc Loebll2009ft . The dif- 
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ferential Thomson optical depth can in principle be de- 
termined from maps of the neutral hydrogen density 
based on redshifted 21cm emission or absorption using 
the same interferometers; very large scale variation will 
be extremely difficult to separate from foreground emis- 
sion. 

7. DISCUSSION AND CONCLUSION 

We have shown that large scale variation in the cos- 
mological baryon fraction is constrained, although not 
strongly, by non-CMB observations. Scatter in observed 
light element abundances constrain cosmological varia- 
tions in the baryon density to be less than 26% (95% 
confidence). Variations are constrained to be less than 
8% (95% confidence) by galaxy cluster measurements, 
provided there are no selection effects at work that lead 
to an artificially low observed scatter in the gas fraction 
compared to the cosmic distribution. 

Within these bounds of variation, there is an opportu- 
nity to ease the tension between the observed Li and the 
standard BBN calculation. In addition, there is evidence 
of an asymmetry in the large scale CMB fluctuations 
which could be caused by differential screening effects. 
However, in both cases the observed effects appear to 
require fluctuations that are much larger than allowed 
by galaxy cluster measurements and well above the up- 
per limits from BBN. It is possible that our Galaxy is in 
an especially baryon-deficient part of the universe, and 
that there is an anomalously large baryon asymmetry on 
scales of the comoving distance to z ~ 7, but this would 



require fine-tuning. 

Variations in the baryon fraction constitute a mech- 
anism for generating B modes in CMB polarization; a 
detection of B modes on large scales would therefore not 
be an unambiguous detection of primordial gravitational 
waves. The existence of such large scale baryon frac- 
tion modulation would be most simply explained through 
an inflation-like mechanism, as there would be super- 
horizon correlations in the large scale baryon fraction. 

In general, cosmological variations in the baryon frac- 
tion (baryon-cdm isocurvature modes) are allowed at a 
level that could be important for galaxy and AGN for- 
mation and evolution, and could also be a foreground for 
future CMB experiments searching for the signature of 
primordial gravitational radiation. Detection of such a 
variation would be extremely interesting for fundamental 
physics, as well. 
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